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ABSTRACT 


fois thesis is an analysis of the one-on-one ASW search 
problem using a random active search strategy in an environ- 
ment that favors the target's counterdetection ability. MThe 
objective is to determine an optimum ping strategy by simu- 
lation of the definite-range problem, approximation by an 


analytical model and use of empirical regression techniques. 
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Tt. INTRODUGEION 


A. THE THESIS Opec TVS 

This thesis documents the analysis of one-on-one ASW 
encounters between a surface searcher using active sonar and 
an evaSive target submarine. The analysis is based on data 
generated by a computer simulation of the relative motion of 
the two adversaries over time. The specific objective of 
this analysis is to prescribe a strategy for selecting a 
searcher ping interval which maximizes the probability of 
detection in an environment which favors the target's counter- 


detection ability. 


B. THE SCENARIO 

A single surface ship 1S assigned to search for a sub- 
marine target of interest using active sonar within a region 
Several humdred thousand square miles in area. The acoustic 
environment is considered homogeneous throughout the area. 
Thus for any particular case, the sonar range is considered 
aw COnStant.. 

1. The Searcher's Tactic 

The searcher's tactic is to move through the area at 

a set speed, changing course randomly at times described by 
an exponential distribution of mean 1/8. The use of the 
exponential distribution for this purpose seems tactically 


prudent because of its memoryless property. The searcher 


pings at times selected at random treomesone distrreurient 
By this tactic, the exact time between succesSive pings is 
made unpredictable. 
ae The Target's Tactic 

The target is patrolling the area of interest at a 
set speed, changing course randomly at times also selected 
from an exponential distribution, not necessarily the same 
as, but independent of, that of the searcher. The target is 
capable of counterdetecting the searcher's transmissions at 
ranges greater than the searcher's detection ranges. It is 
assumed that the target has no method of detecting and 
locating the searcher other than by passive detection of the 
searcher's transmissions. If the target does counterdetect 
the searcher outside the searcher's detection range, he 
Sprints away from the searcher radially at a speed greater 
than that of the searcher. The duration of this sprint is a 
tactical decision made by the target, based on what he con- 
Siders a "safe" range, and can only be estimated by the 


searcher. 


te VObDE ne DESCRIPTION 


It is commonly assumed that the time, Tye required for 
a randomly moving searcher and target to first come within 
some relatively small range of one another is distributed 
exponentially with a mean that is a function of that range, 
the searcher's and target's speeds, and the size of the area 
in which they are confined [Ref. 1]. Intuitively, it would 
seem if the searcher selects a ping strategy that maximizes 
the probability of detection, given that the target is 
within counterdetection range, that strategy tends to mini- 
mize the expected total time, T, spent searching in the area 
for the target. With this in mind, a model of the search 


problem after the realization of To 1s described below. 


a tHe EVENT Disk 

In the model used for simulation, the event disk, l.e., 
me Gegtomewithin am event circle Of radius C, the counter- 
detection range, represents the area in which any interaction 
between the target and searcher must occur. Concentric with 
the event disk is the searcher's detection disk of radius D. 
The event disk is centered on that opponent with the higher 
speed. Figure 2.1 illustrates the case for which the searcher 


is at a higher speed than the target, but the labeling is 


completely arbitrary because of symmetry. If the target is at 





— e 
R: target motion relative tosearcher 


Figumre 2. The Event Disk. 


bee higher specd,~he is placed in the center. All relative 


relationships remain the same. 


Be one EXPECTED SEARCH TIME 
Once the target has entered the event disk, one of three 
Sevemnes MUSE OCCUL: 

(1) The target departs the event disk before the searcher's 
first ping, by virtue of the relative motion between 
the two. 

(2) The target sprints out of the event disk as a result 


of being located in the counterdetection zone but not 
in the searcher's detection zone when the searcher 


pings . 
(3) The target is detected as a result of being located in 
the detection zone when the searcher pings. This event 


completes the search. 
If Events (1) or (2) occur, then there is a possibility that, 
eventually, the target will again enter the event disk. 
Therefore, once the target has entered the event disk for the 
first time, the remainder of the search can be thought of as 
a series of cycles during which the target is either detected 
or not detected. This suggests the use of the geometric dis- 
tribution to describe the process. If the search requires N 
such cycles, then (N-l) of the cycles must have resulted in 
no detection occurring. Therefore, if P is the probability of 
detection, the probability that the search requires n cycles 


Le 


> i oa (Clos) Srrs (2.1) 


Determining, E[T], the expected total time for completing 
the search requires that two additional variables be defined. 


Let 
a (i) —— 2 rs 9! 
Cc 


be a sequence of independent, identically distributed random 
variables describing the cycle time for the i cycle and 


let T, be the time required for the target to enter the detec~ 


d 


tion zone and be detected, given he is located on the perimeter 


Gf the event disk. ween 


sad 
tS See iam a 7 eel 5 (ZED 
i=l 
and 
silane) = ET, ] + (ately pede eM enh ceer 
Removing the condition on n results in the following 
expression: 
EIT} = E{E(T|NJ} = ElT,} + E(t.) } (n-2)(1-p)? 7p 
n=1 
+ E[T,] (2aa)) 
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The summation term is easily collapsed. Let 


ee ea en (1-P)"p 
nel m=O 
= 0+ (1-P)P + 2(1-P)*P +... (2.5) 
Then 
2 3 4 
(l-P)S = (1-P)°P + 2(1-P)~P + 3(1-P) P+... (2256) 


Suebisacting EOuation (226)) from Equation (2.5) yields 


PS = Pp(l1+P)[l1 + (1+P) + (14P)° + ...] ar 


The sum inside the brackets is a geometric series which con- 
verges to 1/P. Therefore 

E(T] = E(T,.] + ({-1)E(T.] + Elta] (2.8) 

0 P © d 
It can be seen that maximizing the probability of detection, 

P, will aid in minimizing the expected total search time. It 
is for this reason that this thesis concentrates on the prob- 
lem within the event disk; that is, finding a ping rate which 


Maximizes P. 


el: 


C. CRITICAL ASSUMPTICNS 
Several assumptions have been made to simplify the model 
and analysis of the generated data. In addition to those stated 
previously, the following also apply: 
(1) The occurrence of detection and counterdetection events 
1s determined using a definite-range or "cookie 
cutter" model. 


(2) The target and searcher have negligible barfle areas. 


(3) Detection and counterdetection ranges are not degraded 
with speed. 


(4) There is no convergence zone considered, nor are there 
any GapS In sehesecVentwenmerc, 


(5) The target is strictly evasive. 

(6) Counterdetection range as used in the model should be 
considered the target's minimum desired range to the 
searcher. It is assumed this range is at least twice 
the detection range. 

D. THE REQUIREMENT FOR A PING STRATEGY 
The assumed existence of a definite counterdetection range 


greater than the searcher's detection range requires that the 


searcher have a well-defined minimum interval between any two 





active pings. This minimum interval is: 
» Cc —op 
Tnin ~ U+V ae 


where C is the counterdetection range, D the detection range, 
and the denominator is the sum of the two speeds. This is 
merely the time required for the target to move from the 
perimeter of the event disk to the detection zone at the maxi- 


mum attainable relative speed. 


ie? 


There is also a maximum practical ping interval that is 
not as well defined. That is, if the searcher pings very 
infrequently he loses the opportunity for detection because 
the eaves may transit in and then out of the detection zone 
between pings. This "maximum" should depend upon the size of 
the detection zone, and the relative motion between the 
searcher and target. It is because of the relative motion 
aspects that this maximum practical interval cannot be defined 
as easily as the minimum. The existence of a minimum ping 
interval below which the probability of detection is zero, 
and a "maximum" ping interval beyond which the probability of 
detection is small, implies that the probability of detection 
reaches a maximum between these two extremes. This maximum 
should be a function of the ranges and speeds specific to each 
particular case. Therefore the first step is to determine a 
ping strategy as it depends upon the independent variables. 
This will be done by simulation, approximate analytical modeling, 


and a blending of the two by an empirical regression technique. 


E. DATA GENERATION 

The equation for probability of detection if both searcher 
and target remain on constant courses and speeds is complex 
but can be solved using polar coordinates and some trigonometry. 
However, in the problem at hand, both relative speed and its 
direction change randomly. This urges the use of simulation 


to generate data on relative courses and positions. 


i¥s 


The simulation program for the definite range problem 
(included in Appendix C for informational purposes) supplies 
as output the number of detections, the number of counterde- 
tections, and the number of times the target departs the event 
disk before the first ping, for a given ping interval, detec- 
tion range, counterdetection range, searcher speed, and target 
speed. The frequency of course changes, the searcher course, 
and the target course are determined by random number generators. 
All relative motion is placed on the target using trigonometric 
arguments with the searcher remaining at the coordinate origin. 

Because each interaction begins with the range between 
the searcher and target decreasing to C, it follows that the 
initial relative velocity must be directed into the event 
disk. To accomplish this in the simulation, the first rela- 
tive velocity vector was determined by assuming that the 
searcher's speed component directly toward the target was just 
greater than the target's speed. After the initial leg, all 
motion was unconstrained. So after several course changes, 
the effect of the initial leg is "forgotten" by the process. 

An alternate method, not used in this thesis, would be to 


let the target's initial relative angle on the bow be selected 


fremaa. cosinesd SEG be roOne Taito 
TL 
P (Be waiewn oy AO < eam cos 6 dé 
—-17/2 
= F(sin 9 +1), -1/2 < 6 < 1/2 (2.10) 
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This procedure is suggested by Koopman's observation that 
when searching for a stationary target, the bearing of 
initial detection will have a cosine distribution [Ref. 1]. 

It is not known which of these methods is more correct. 
It may be possible to derive an exact expression for the joint 
distribution of the target's relative speed and course at 
detection, but this was not accomplished here. 

Once the time for the first ping is reached, the replica- 
tion is stopped. The range between the two is calculated 
and compared to the values for C and D to determine which of 
the possible events has occurred. The outcome is stored and 
the whole process is repeated for the desired number of trials. 
The data is then analyzed to obtain estimates of the proba- 
bility of detection and other relevant quantities, such as 


expected times within and without the event disk. 


il 


III. DETERMINATION OF THE REGRESSION MODEL 


Once the data has been generated, regression may be used 
to determine an empirical predictive formula for the optimum 
ping rate. Armed with only intuitive hypotheses, the search 
for the "best" functional form of the input variables would 
be difficult. So it is desirable to find some theoretical 
guidelines for selecting candidate explanatory variables to 


use in the regression model. 


A. USE OF THE VON NEUMANN FUNCTION 

The theoretical model selected for use is one for energy 
transmission and return when the target's motion is modeled 
by a diffusion process. Define the searcher's location as 
the origin on a Cartesian plane, and define the target's 
location at time t, as [X(t),Y(t)]. Let the target undergo 
Brownian motion so that its location at time t is described 


by 
Z 
XC) SNe s(0 oe 
(3 cy 
2 
VCC) See (Y (On, © ep) 
and let the searcher's probability of detection, for a ping at 


time t be 


eG 


Plt) Meee (t)) = e4 


where R(t) is the range from the searcher to the target and 


ac is a, thus far, unspecified constant. The constant 6 in 
the detection function (Equation (3.2)) plays a role analogous 
to that of the detection disk radius, D. In particular, 6 is 


1/2 . 9.607. 


Brae eeange where the probabulity ot detection is e 
Removing the condition on position by integrating over 


all values of X(t) and Y(t) results in the following expression 





Ota (t): 
© =S(K(t)“)/6o | -F 1 ((X(t)-X(0)) 7) /o*t] 
P(t) = f fe e ax 
_ donee 
“S((¥(t)“)/6° “SL ((¥ (t)-¥ (0) ) 7) /o*] 
e n= dy 
V2nO-e 
-F1(R(0)°/87) / (1+0°t/6*) ] ; 


=> i ——— 
UW AGE sae) 


Usangumence mactural log function to linearize P(t), and 


differentiating with respect to t yields 


Dp Z 2 
ol 0es Ge a 0/6 ee eo Eh (3.4) 
= a) Ome. 2 , 
l1+o t/6 6 . 


Setting the derivative equal to 0 results in the following 


expression for T*, the optimum ping interval: 














iL 
i (/R(0) ; 1 
Z 52 
mx = 5 (3350) 
(2) 
6 
Therefore 
2 
i =O 1 T = ey ek 
> = 
6 
and 
ie 2 
> eo) - l 
Ta aa otherwise (Soy 
(=) 

Referring to Equation (3.2), suppose energy from a ping 
transmitted from location (0,0) reaches the target with 
PLObaba lity 

RAGE). 2 
a5 (=e) 
O 
e ; (30 


and the probability of the signal returning to the searcher 
given it reached the target is 
3 ie eZ 


il 
=> ae 
e 7 (3.8) 
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Then the probability of the searcher detecting the target 


becomes 





ey Ree ob R(t), ¢ 
Paes DS 
P(detection R(t)) = e : e it 
a 
~5R (t) ee 
O aL 
= e 
-1 R(t), 2 
= ne (3.9) 


which is also a Von Neumann detection function. 

The intention at this point is to use the form of Equa- 
tion (3.5) to develop an approximation for the optimal ping 
interval TOR for the definite range problem. This is accon- 
plished by setting R(0) in Equation (3.5) equal to the 
counterdetection range, C, and recognizing that the parameter 


§6 is a “characteristic range" for the Von Neumann function. 


SO 


+ = eS) = Bess (3.10) 
6 6 6. R is 
O a 


where R and r are the ranges associated with the outbound 
and inbound acoustic paths. 


SUDSterPeuemnGg InitoeEGuation (3.5) yields 


teeth, +) - ji 
yA ce re 
* — pein B= 
TDR 7 ors il (pmo 
o (— + —) 
2 2 
ei Yr 


Letting r now be the active detection range, DD; restmaeosin 


the following expression: 


i Cac" sy pee 2/C* 
TOR oe) x1 2 2 
20 1/C, al Db 
| Gee 1 p- =a ee 


If 


] 
oe) 1/D- eee 


| e= ce | 
20 c*/D + 1 


] (3). eZ) 


B. THE EXPLANATORY VARIABLES 

It remains to replace the diffusion constant o with appro- 
priate random tour model parameters u, v, and X. Lambda is 
the parameter for the exponential distribution describing the 
minimum time between course changes for either of the adver- 
Saries. If 67 is a diffusion rate describing the relative 


motion Of the ewe, Enenm 
5 a (2mies 


is dimensionally correct and has some theoretical appeal. 
Specifically, an unconstrained two dimensional random tour 
with constant speed V and rate of course change Aa can be 
approximated for large times by a diffusion process with 


constant v°/x, (Ref. 2]. When two particles are simultaneously 
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conducting random tours, then the composite random tour in 


relative space has a rate of course change: 
= (se 4) 


and a random speed Sp: If the angle between the V and U 
velocity vectors is uniformly distributed between 0 and 27 


radians then the expected value of sé 1S: 


27 
= i (U? +v* -2 UV cosé)dae 
0 


2 
E[S,] 
=e Pe ema (or) 


Using (y? ARGS) aS a representative squared speed for the 
composite random tour and A for the rate of course changes 
maclads Equatien (3.13). 

Fguation (3.12) is related to, but probably unequal to, 
the optimum ping interval for the definite range law model 
that is being simulated. In order to better adapt the Von 
Neumann model results to the definite range law data, one 
can redefine TAR as 
ee 7-3) 
Ge We) ae feo + 1 


ze 
T*, = B (ene 


DR i! 


and determine the parameters Bye B and B, by regression, 


Das 


using as explanatory variables, the quantities 


zal 


ec Coo = sl 
lr and | 
(U~ +V") C/Dawe*1 


ts is acknowledged that Equation (3.16) is probably not 
the "best" definite range law extrapolation of the Von Neumann 
result of Equation (98)). In pamewewlar, sectemmne Os equal 
to the detection range D seems suspect since the active 
detection process involves two-way propagation and Fi considers 
primarily the return path. Nonetheless, Equation (3.16) was 
tested as a candidate regression model and, as the next section 


documents, it performed quite well. 
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iVo THESRESUETS 


A total of 180 cases were run uSing 4 different speed 
combinations, detection ranges from 3 to 16 miles, and several 
counterdetection ranges between 20 and 40 miles. The proba- 
bility of detection was measured at each of 50 different ping 
intervals, beginning at the practical minimum and stepped in 
0.05 hour increments. A total of 500 trials were run at each 
ping interval. The observed T* was that ping interval at 
which the maximum probability of detection occurred for each 
case. In the event of a tie, the earlier time was used. 
Figure 4.1 illustrates the dependence of the probabilities of 
detection, counterdetection, and departure on the value of 
ping interval. As expected, the probability of detection is 
O until the ping interval is greater than (C-D)/(U+V), in- 
creases to a maximum and slowly decreases to a small positive 
value. The probability of departure increases with ping 


interval, while the probability of counterdetection decreases. 


ye  fHE BMPLRICAL PREDICTION FORMULA 
Performing linear regression on the data, using the 
explanatory variables discussed previously, produced the 


following equation for predicting TAR? 


Z Cpe 2 
XC :0-661 °C /D LO fey 3 (am 


yO) ] 


A = 0.74[ 7, 
(U~ +V°) 6. /D +41 
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This formula explained 89% of the total variation in the data 
and all coefficients were significant at the 0.99 level using 
Student's t statistic. 

Tables (1)-(4) of Appendix A list the input variables, 
observed T* and predicted T* for each case. It is immediately 
apparent that the probability of detection at the predicted 
value of T* is consistently less than or equal to that at 
the observed value of T*. This is because of the initial 
relatively crude method used to choose the observed T* for 
each case. The values for the difference in the probabilities, 
a measure of the prediction validity, are relatively small. 
Figure 4.2 is a scatter plot of the probability of detection 
at the predicted T* versus that at the observed T*. The ideal 
plot would be a straight line of slope 1 through the origin. 
The least squares fitted line througdithe data has a slope of 
0.94 and intercept of -0.01. Many of the points away from the 
Pagetaumeanebe Cxplarmed by examination Of the raw data. 
Often, the predicted and observed optimum ping intervals are 
within 0.10 hours of one another, but the variance of the 
sampled binomial distribution causes the two detection proba- 
bilities to appear farther apart than might actually be the 
case. It should be mentioned that smoothing half of the raw 
data using running medians followed by running averages 
(Hanning) and using the same regression model did not alter 
BiemCOoEtTetomis OF the predretion Lormula significantly. 


It did, as expected, tend to decrease the difference between 
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the value for probability of detection at the predicted T* 


and that at the observed T*. 


Bee ecoNCLUSIONS 

Rie Most signuammecanterinding @seehat it is possible to 
determine an optimum ping strategy within the limitations of 
the model. If a searcher were tasked with such a search, he 
would not ever want to ping more frequently than the minimum 
ping interval (C-D)/(U+V). This assumes the searcher has a 
good estimate of the counterdetection range and speed of the 
target. On the average, he would want to ping at a rate that 
is slightly less than that prescribed by the prediction. 

Where arceetwestreasons@:or this: (1) the probability of 
detection decreases more slowly to the right of T* than it 
iicreases Atwemeomlerlt- “ange(2)many deviation to wehe right of 
T* merely increases the probability of the departure event 
occurring. The departure event is considered less detrimental 
to the search effort than the counterdetection event, for 
which the probability of occurrence steadily decreases. 

The Von Neumann function seems well-suited as a theoreti- 
cal foundation for determination of an optimum ping strategy. 
This suggests that research into its use in a more realistic 
model of active sonar search might prove valuable in predictions 
Cf tls somes ine Cxtremely sharp cut offpyof the) definite 


range law is certainly an artificiality. 
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APPENDIX A 


TABLES 
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COMPARISON OF OBSERVED AND PREDICTED a 


Se 


KTS (TARGET SPEED. 
LAMBDA: 1.5/HR 


5 


SEARCHER SPEED: 


F MOPAR OO MORNE We VOLDOOOOOODO Monasanoowaon 
A* NUHNNHHHOW 


fx} 
WV) 

fQ SOP IOA RA BOE SOB OO0ONONDOO0OiNINOCOCCOuMMMOn 
OF MOMMENOAANNSROAMER 


o~—, 

= 

ae PMNS DMNOIONMNO~ DN HNO OHO MOM CIM HOMOricieo dune 
5 rd odode mood 
ane, 

= 

Ze OC N Lf) ™ oO O 

oy N N N N ~ sH 


28 


VERGE S 
eel S 


TABLE 2 
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TABLE 3 
COMPARISON OF OBSERVED AND PREDICTED T* 
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TABLE 4 
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